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Crossovers Between Elastic, Plastic and Quantum Creep in 2D 
YBa 2 Cu307/PrBa 2 Cu 3 07 Superlattices 

X.G. QiuO, V.V. Moshchalkov, and Y. Bruynseraede 
Laboratorium voor Vaste-Stoffysica en Magnetisme, Katholieke Universiteit Leuven, Celestijnenlaan 200D, B-3001 Leuven, 
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G. Jakob and H. Adrian 
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Two-dimensional (2D) vortex dynamics was studied in an YBa2Cua07/ PrBa2Cu3C>7 superlattice 
by measuring the I-V characteristics. In the high current limit, 2D collective creep was observed with 
an activation energy characterized by U(J) oc j _M . A dislocation mediated vortex melting happened 
when the temperature increased. In the low current limit, the exponential growth of energy barrier 
for elastic motion was prohibited by the plastic deformation of vortices. A plateau in the resistive 
transition was observed, which was attributed to possible quantum tunneling of vortices. Our results 
suggest that a 2D vortex glass can not exist at any temperature, including T=0 K. 
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The unusual resistive and magnetic behavior of high 
T c superconductors (HTSCs) has greatly stimulated the 
efforts to understand the impact of disorder, thermal fluc- 
tuations and dimensionality on the vortex dynamics [0. 
The theory of collective flux creep describes the thermally 
assisted motion of vortices in a random potential caused 
by quenched disorder . Central to the theory of collec- 
tive flux creep is the elasticity of the vortex system. Due 
to the finite values of compress modulus en and shear 
modulus C66, when moving from one metastable position 
to another one, the vortices have the tendency to jump in 
bundles with a radius R c in order to balance the elastic 
energy and the energy related to Lorentz force. The en- 
ergy barrier for activation U c increases exponentially with 
a decreasing current because R c gets larger and larger. 
However, this picture breaks down when plastic motion 
of vortices is taken into account . For a 2D vortex lat- 
tice, due to the finite energy for creation of dislocations 
or dislocation pairs, a plastic motion of vortices is favor- 
able when the energy barrier exceeds a certain threshold. 
The plastic creep sets a cutoff for the exponential growth 
of U c . The plastic behavior of vortex system has been 
studied by several groups Q . Recently, plastic motion of 
vortices was directly observed on Nb films with artificially 
introduced pinning centers by using Lorentz microscopy 

Quantum tunneling of vortices is expected to replace 
thermal activation as the dominant dissipation mecha- 
nism when the temperature is low enough . Tunneling 
of vortices has been observed in transport and magnetiza- 
tion measurements fip|,p| Ji^JTl[ | . It was found that the re- 
laxation rate of magnetization did not extrapolate to zero 
when the temperature went to zero, and this phenomenon 
was attributed to quantum tunneling of vortices |Q . Vari- 
able range hopping of vortices was observed on ultrathin 
Pb films H and YBCO single crystals with columnar de- 
fects Ephron ct al. |ll]] observed a transition from 



thermally activated to a temperature independent resis- 
tance on Mo4 3 Ge 5 7 ultrathin films, and this resistance 
was attributed to the quantum creep of vortices. 

Plasticity of vortex lattice (VL) and tunneling of vor- 
tices are two possible contributions prohibiting the exis- 
tence of vortex glass state which is important for future 
applications. Up to now, they are not well studied and 
our understanding to these contributions is still limited. 
This strongly motivates further studies of the interplay 
between elasticity, plasticity and tunneling of vortices in 
a 2D vortex system by utilizing suitable samples. In this 
paper, we report our study of the vortex behavior in 2D 
YBa 2 Cu 3 07/PrBa 2 Cu 3 07 (YBCO/PBCO) superlattice 
at different temperature, magnetic field and excitation 
current. We observed a crossover from collective (elastic) 
creep to plastic creep of vortices when the elastic limit 
was reached for a VL, which means a 2D vortex glass can 
not exist at finite temperature. A dislocation mediated 
melting of VL was observed in the high current limit. A 
possibility of quantum creep of vortices at low temper- 
ature and thus the non-existence of 2D vortex glass at 
T=0 is discussed. Our results provide for the first time 
a comprehensive picture for the vortex dynamics within 
a 2D VL in artificial superlattices. 

The c-axis oriented [YBC0(2AA)/PBC0(IUA)} 25 
superlattice was fabricated by in-situ magnetron sputter- 
ing jL2j. The X-ray diffraction showed satellite peak up 
to the third order, indicating the high quality of the sam- 
ple. This film was photolithographically patterned into 
strip with a width of 0.1 mm. DC I-V characteristics 
were recorded with currents generated by Keithley 238 
and voltage detected by Keithley 182. Resistive transi- 
tion /? ac (T) was measured by four-terminal ac locking-in 
technique with an excitation current of 1 /xA at a fre- 
quency of 17 Hz. The magnetic field was generated by a 
15T Oxford superconducting magnet. During the mea- 
surements, the field was kept perpendicular to the film 
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plane and the temperature stability was better than 10 
mK. 

Shown in Fig. 1 are the representative I-V characteris- 
tics measured at 2 T with the temperature varied from 2 
to 32 K. The I-V curves were taken with a temperature 
interval of 0.1 K. Each I-V curve was obtained by aver- 
aging over four individual curves with current applied in 
positive and negative directions. Two features are clearly 
visible: (i) in the high current limit (j > 5 x 10 3 A/cm 2 ), 
I-V curves show downward curvature at low tempera- 
tures, changing to upward curvature 
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FIG. 1. Selected I-V curves at different temperatures rang- 
ing from 2 to 32 K with an interval of 2 K. Inset: temperature 
dependence of p at 11=0.5 T and 2 T, respectively. 

at high temperatures; (ii) in the low current limit, linear 
I-V curves were observed for all the temperature region 
over several order of magnitude for current and voltage. 
The linear resistivity ph extracted from the I-V curves 
are shown in Fig. 2 together with the R-T curves p ac (T) 
in the Arrhenius plot. We find that pl coincides with 
p ac quite well. In the high temperature part, a thermally 
activated flux flow (TAFF) behavior is clearly visible, 
when the temperature decreases (<10 K), the resistivity 
deviates from the TAFF behavior, and reaches a constant 
value in the lower temperature at each field. 

Due to the large PBCO layer thickness, YBCO layers 
are essentially decoupled, so each YBCO layer can be 
treated independently. Moreover, each YBCO layer con- 
sists of only two unit cells (24 A), which is much smaller 
than L c (— lOOA), the coherence length of vortices along 
the c-direction. So we treat the whole vortex system as 
2D. To understand the above I-V characteristics in the 
whole temperature and current ranges, we discuss the 
following 2D vortex dynamics scenarios. 

2D collective creep and vortex melting- It is expected 
that in a 2D disordered superconductor, 2D collective 
creep of vortices plays an important role in the dissipa- 
tion process. 2D collective ceep theory is applicable when 
R c = a (eod/U p )(£/2a ) 2 > a w ($ /S) 1/2 , the vortex 
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FIG. 2. Temperature dependence of R ac at difference ap- 
plied magnetic field. The solid lines are experimental data. 
From left to right: H=0, 2, 4, 6, 8 T. Fitted circles are 
data points of Rl extracted from the linear parts in Fig.l. 
Inset: Activation energies from a TAFF fit to R ac (T), 
filled circles are the obtained data, solid line is a fitting by 
(7=141. 8-105.9/nB. 

spacing, where eq = & 2 /16ir 2 \ 2 is the energy for a pan- 
cake vortex per unit length, U p the pinning potential, d 
the length of the vortex along the field direction, £ the 
coherence length |^3| . From Fig. 2 we have U p « 100/V 
with d=24 A, A « 1400 A, f » 15 A, a Q «500 A, and 
R c -5000 A w 10a . Therefore we think that 2D collec- 
tive creep theory is suitable in the present study. When 
vortices undergo collective creep, the activation energy 
is U(j) oc (j/io)~ M ) where jo is a characteristic current 
density. Within the TAFF approximation, the I-V curves 
caused by the presence of the energy barrier U(j) is [|| 

E = E e-^V/ jo ^ . (1) 

According to Eq. (1), the logarithmic plots of the I-V 
curves with downward curvature were fitted by a power 
law, with IhEq, (Uo/kT)jQ , and p as fitting parameters. 
The obtained exponents p are shown in the inset in Fig. 
1. The typical values for lnE and (Uo/kT)jQ are —10 
and —1000, respectively. The p, values change gradually 
from 0.8 to 0.2 as the temperature increases. This result 
agrees well with that of Dekker et al. who found the same 
trend in the varoation of p upon changing the tempera- 
ture Jl^l . We interpret the change in values of p as fol- 
lows: as pointed by Vinokur et al. jl3[| , depending on the 
current density, temperature and magnetic field, the bun- 
dle size will be different and consequently p will take var- 
ious values. For vortex bundle with medium size (defined 
by jib w ]c {Rl/Ka 2 f/^ < j < Jmb » ]c (R c /A) s / 15 ), an 
exponent /i=13/16« 0.8 is expected. When j < jit, the 
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bundle size involved in a jump increases, and the expo- 
nent p changes to 1/2. We notice that when the magnetic 
field is fixed, the bundle size is R c cx £, 2 /\ 2 U p . At tem- 
peratures not too close to T c , where £ and A are nearly 
constants, U p decreases with increasing temperature, and 
therefore R c will increase. As a result, the bundle size 
will change from medium to large thus decreasing p value. 

As the temperature increases, dislocation mediated 
vortex melting occurs when the following condition is sat- 
isfied: m 



Ac§§a\djkT — Air 



(2) 



where A ~ 0.4-0.7 is a parameter due to the renormaliza- 
tion of C66 from the defects in the VL and the nonlinear 
lattice vibration. We identify the melting temperature 
as the one at which the I-V curves at large current show 
a power law dependence, i.e., V oc I 13 . For H=2 T, we 
obtain V oc I 3 at T=12.87 K. The shear modulus cqq has 
been calculated by Brandt Jlq] who showed that 



c 66 = [B c2 (t)/4p ]b(l - 0.296)(1 - bf 



(3) 



where B c (t)=B c (0)(Ti 2 ) is the thermodynamic critical 
field, t=T/T c and b=B/B C 2(t). Inserting the typical val- 
ues of B c and B c2 for YBCO with B c (0) wlT and B c2 (0) 
wlOO T, together with B =2 T, T=12.87 K, and T c =31.5 
K, we obtain A w 0.45, which agrees well with the pre- 
viously obtained value Q| . 

Above the melting temperature, the vortices are in a 
liquid state pinned by residual forces. The I-V characeris- 
tics can be well described by the classic Anderson-Kim 
model with E — Eq sinh(j/jo)- 

We notice that Dekker et al. [[l7j did similar work on 
ultrathin YBCO films. In their study, only I-V curves 
with upward curvatures were observed and the glass cor- 
relation length £vg diverges at T=0. They concluded 
that a 2D vortex glass could not exist at any tempera- 
ture T>0. In our case, we did observe a crossover from 
downward to upward curvatures. This discrepancy could 
be possibly due to different length scale involved in the 
experiments. As we know, when applying a current, a 
vortex pair of a characteristic size Lj = {ckT / j^o) 1 / 2 is 
excited |fl8f . The current probes the vortex dynamics in 
a length scale of Lj. With j ~ 5 x 10 3 A/ cm 2 , T~ 10 K, 
we get Lj ~ 1000 A < R c . Therefore, we were probing 
the collective creep behavior in the high current limit. 
While in their case, the length scale was probably larger 
than i? c , and thus vortex glass behavior was detected. 

Plastic motion in the low current limit- The collective 
pinning theory depends critically on the elasticity of vor- 
tex system. Due to the finite energy for the generation 
of dislocations or dislocation pairs in a 2D VL, an infi- 
nite energy barrier in the elastic limit is not expected. 
When the energy barrier reaches the characteristic en- 
ergy for a small dislocation pair, it will be cut off and a 
crossover from elastic to plastic motion should occur || . 



The motion of dislocation pairs can be well described 
by TAFF with an activation energy nearly independent 
of the current density. Therefore, linear I-V curves are 
expected in the plastic region. As a result, there is no 
vanishing linear resistivity, and thus a 2D vortex glass 
does not exist at any finite temperature. This is exactly 
what we observed in our I-V curves. We see that at a 
current density about j~5x 10 3 A/cm 2 , the I-V curves 
with downward curvature change to linear at low current 
densities. This linearity is valid in the range of more than 
two orders of magnitude in the current and voltage. The 
typical energy for a dislocation pair made up of two edge 
dislocations separated by qq is [0 



16tt 3 A 2 (T) 



HBo/B) 



(4) 



With A «1400 A, we have U e « 500 K. This is consis- 
tent with what we have obtained from a fitting to the I-V 
curves we made earlier. With {Uo/kT)j^ ~ 1000, jo ~ 
10 5 A/cm 2 , we have Uq ~ 300 K, which agrees quite 
well with U e . This explanation is further supported by 
the InB dependence of activation energies which were ex- 
tracted from the slopes of the Arrhenius plot as shown in 
Fig. 2. So we conclude that in the low current limit, the 
dissipation is governed by plastic motion of dislocation 
pairs in the VL. 

The above picture is consistent with Monte Carlo sim- 
ulations of the VL subject to random disorders under a 
driving force p9| . It is found that at large driving force, 
the vortex motion is elastic. As the driving forcing force 
decreases, the VL becomes defective, and the dissipation 
is dominated by plastic motion. 

Possible quantum tunneling of vortices at low T- as 
we can see from Fig. 2, below T~10 K, the resistivity 
gradually deviates from the activated behavior, showing 
a level off toward constant values (^0|. The resistivity 
at the plateau grows with the increasing magnetic field. 
To get such a plateau within the framework of TAFF 
would require an activation energy Uq which increases 
linearly with T, i.e., Uq oc oT. We are not aware of any 
mechanism that can provide such kind of an activation 
energy. Same resistive plateau in the R-T curves was also 
observed by Ephron et al. on Mo43Ge57 thin films. 

Generally, it is believed that quantum tunneling, if it 
does exist, dominates the dissipation at temperatures 
below several Kelvin and that the resistivity from a 
quantum creep is nearly independent upon temperature 
HUmB We are not sure why we observed quantum 
tunneling behavior at such a high temperature where the 
thermal energy of vortices is quite high. One possibil- 
ity is a large normal resistivity p n and a small coherence 
length £. The tunneling rate 7 is determined by the ef- 
fective Euclidean action for the tunneling process, 



7 oc exp(— S e E f /h), and S E /Ti in turn is proportional 
to (fi/e 2 )^/ p n pl[ . Therefore, tunneling is favored by a 
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small coherence length £ and a large normal state resis- 
tivity p n . Experimentally, Ephron et al. reported that 
they were only able to observe such a plateau on thin 
films with large sheet resistance, this is consistent with 
theoretical consideration. 

The existing results on the temperature To at which 
the dissipation shows a crossover from thermally acti- 
vated to a quantum tunneling range from hundreds mK 
to tens Kelvin. Ivlev et al. |2^] showed that, for YBCO 
compound, a crossover from thermally activated to quan- 
tum tunneling dissipation happens at To ~ 50{j c — 
j) 1/2 /jc /2 K. Meanwhile, Stephen || found that T «2 
K. Experimentally, a crossover temperature ranging from 
3.5 K to 10 K was reported for TIBaCaCuO single crys- 
tals ||. When vortices are localized by columnar defects, 
a crossover temperature from collective creep to quantum 
creep (variable range hopping) as high as 35 K was re- 
ported in the high current limit (^10 6 A/cm 2 ) Jicf| . It is 
worth noting that Blatter et al. (2l]] found the suppres- 
sion of the tunneling rate for increasing field. In this case 
the interaction between voritces increases, the viscocity 
grows and more vortices are participating in a single tun- 
neling event. This is consistent with our observations. 
We found that the temperature, at which the resistivity 
starts to deviate from TAFF behavior, decreases with the 
increasing magnetic field. 



or an interstitial in a VL, and a 2D vortex glass can not 
exist even at T—0 K. At present, due to the lack of a 
suitable theoretical model, we are not able to do a quan- 
titative analysis of the appearance of the puzzling plateau 
in the low T region. Further theoretical and experimental 
work at that subject is greatly encouraged. 

The main conclusions of this work are summarized in 
fig. 3. We found that the vortex dynamics is essentially 
length scale dependent. In the high current limit, 2D col- 
lective creep was observed which is elastic in nature and 
characterized by an energy barrier with an exponential 
growth with current density. The 2D VL undergoes a 
dislocation mediated melting at a temperature T m . In 
the low current limit, the exponential growth of energy 
barrier is cut off by the generation of dislocation pairs in 
the VL and their subsequent motion or quantum tunnel- 
ing. The linear resistivity at low temperatures (<10 K) 
is possibly governed by quantum tunneling of pancake 
vortices. With increasing Lorentz force, the driven VL 
recovers its elasticity at high current density. 

This research has been supported by the Belgian IUAP 
and Flemish GOA, FWO and Bilateral Flemish-Chinese 
(BIL 97/35) programs. 
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FIG. 3. A schematic j-T phase diagram for a 2D VL at 
a certain magnetic field. The line T m is determined by Eq. 
(2), and line T p by Eq. (4). In the quantum creep and plas- 
tic creep regions, the linear resistivity does not vanish. The 
collective creep and Anderson-Kim creep regions are charac- 
terized by I-V curves with E oc e 3 ^ and E oc sinh(j/jo), 
respectively. 

Therefore, our results suggest the possibility of tunnel- 
ing of dislocation pairs which are eqivalent to a vacancy 
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